Testing whether muon neutrino flavor mixing is maximal 
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The small difference between the survival probabilities of muon neutrino and antineutrino beams, 
traveling through earth matter in a long baseline experiment such as MINOS, is shown to be an 
important measure of any possible deviation from maximality in the flavor mixing of those states. 
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How really close to maximal Q is the flavor mixing 
of muon neutrinos and antineutrinos discovered in at- 
mospheric neutrino experiments We propose in this 
Letter a way of probing the deviation, if any, from this 
maximality. The idea is to measure the small difference 
AP W = P aa - P m with P MM ee P [^(0) -» v M (L)] and 
Ppp, ee P [Pji(0) — * P/j,(L)] representing respective sur- 
vival probabilities of muonic neutrino and antineutrino 
beams after passage through a distance L in earth mat- 
ter. Suppose we neglect the subdominant oscillations 
due to the smaller mass difference relevant to the so- 
lar neutrino problem, but retain a nonzero mixing be- 
tween the first and the third family neutrinos. Working 
to the linear perturbative order in the earth matter ef- 
fect, we are then able to show that AP MM is proportional 
to \U e3 \ 2 (l - 2 1 [/^ | 2 ), U being the Pontecorvo-Maki- 
Nakagawa-Sakata neutrino flavor mixing matrix, with a 
computable proportionality coefficient. Thus the effect is 
linear in the deviation l/y/2— |£7^3| from maximal mix- 
ing. 

Much has been learnt recently about neutrino masses 
and mixing angles from solar js| , atmospheric 0] , re- 
actor 0, U and long baseline [a] studies [3] . We know 
now that the squared mass difference between one pair 
of neutrinos comprising and v r , mixed nearly maxi- 
mally ((9 23 45°) @ is |<5m| 2 | ~ 2 x 1CT 3 eV 2 . We also 
know that the squared mass difference between another 
pair of neutrinos, involving v e mixed by a large angle 
(#i2 ~ 30°) with a nearly equal combination of and 
v T , is <5to 2 i ~ 7x 1CP 5 eV 2 . Furthermore, the mixing be- 
tween the third possible pair, characterized by the angle 
#13, is known to be quite restricted, as elaborated below. 

Yet longer baseline experiments with v, v 
beams/superbeams promise to be the wave of the 
future in neutrino physics. MINOS [|J and off-axis 
NUMI 0| are forthcoming experiments and will be 
followed by JPARC [ll|, CNGS G3 and other efforts. 
Many theoretical and phenomenological analyses 13] of 



physics issues pertaining to long baseline experiments 
have been carried out meanwhile. However, the focus 
of recent studies has largely been on appearance experi- 
ments: specifically, the 'golden' (y e — > u^) and 'silver' 
(vp — > v T ) channels and their synergy in probing leptonic 
CP violation as well as the mixing between the first and 
third family neutrinos. Less attention has been paid to 
survival probabilities, the originally measured quantities 
|l4| in atmospheric neutrino studies. This is since 
they do not yield direct information on those aspects. 
However, it is the survival probabilities for and 
beams/superbeams, specifically their difference, which 
should help determine the deviation from maximality, if 
any, of the flavor mixing of muon neutrinos. 

Earth matter directly affects only neutrinos with clcc- 
tronlike flavor. However, this gets induced into the other 
neutrino flavors indirectly through mixing cum oscilla- 
tion effects. For neutrinos of muonic flavor, there are 
two sources of such an occurrence: (1) mixing between 
the electronlike flavor and the third physical neutrino 
through the factor \U e3 \ 2 and (2) subdominant oscilla- 
tions between electronlike and muonlike neutrino flavors, 
driven by Sm^i- It is generally known, and we will show 
later, that the effect of (2) is quite small for the baseline 
length of MINOS . This means that one can take the 
earth matter effect in such an experiment to be due to 
(1) only. We also assume that the actual value of |C/ e 3| 2 
is not as small as one order of magnitude less than its 
current upper bound 0.05 — 0.07 [15j so that our effect 
will be measurable. Moreover, we shall treat this effect 
to the lowest perturbative order 0] in A = \/2GfN g , 
N e being the average electron density of the earth. This 
will also be justified later. On the other hand, if |£/ e 3| 2 
turns out to be smaller by more than an order of mag- 
nitude than its current upper bound, the content of our 
Letter will prove empty. 

The effective Hamiltonian for neutrino oscillation in 
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matter can be written as 
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Here = (m| — m 2 )/ (2P„), rrii and P^ being the mass 
of the «th (physical) neutrino and the beam energy re- 
spectively. For simplicity, let us work in the uniform 



earth density approximation, though our results extend 
to the more general variable density case, as shown by 
using the evolution operator formalism [l7j. When A21 
is neglected in comparison with A31, the effects of the 
solar neutrino mixing angle 612 and of the CP violating 
phase 8 in U become inconsequential. Then, in the stan- 
dard parametrization and with cy = cos 9ij , Sij = sin 9ij , 
H takes the form 



H ~ 
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The RHS of @ can be diagonalized by the unitary matrix 
U, yielding the eigenvalues Ai.,2,3, given by 

h = 0, (3a) 

«i,3 = ^(A 31 +A T B), (3b) 

where 

B = VA31 2 + A 2 - 2A 31 Acos29 13 . (4) 
Moreover, U can be written as 0] 
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where the angle 9 m is related to #13 and A by, 

A31 sin 2#i3 

tan 29 m = ~r ^ 7 

A31 cos 26*13 — A 

The survival probability, after propagation through 
a distance L in matter, is [18j 

P w = l - 4(|C/ All | 2 |^ 2 | 2 sin 2 ^-iiL 



2|tt |2 -2 ^ - ^1 



+ |t/ M irl^3rsm 



L 



+ |^2| 2 |^3| 2 S in 2 ^i 



1 - -i I S23 C 23 S L sin 2 + f —L 



s 2 3 s e m c 9 m sin 2 



BL 



222 -2 A31 + A + B , 
c 23 s 23 c e m sm i . (7) 



Ppp can be obtained from P w simply by changing A to 
—A 0| . Moreover, in vacuum with A = and hence 
?i = = I2 and l 3 = A31 as well as U = U, we have 



= P m = 1 ^ 4 I^3| 2 (1 " l^3| 2 ) sin 



A 31 L 



(8) 



The oscillation probability for a v^jv^ traveling in vac- 
uum, namely 1 — P^ c , is maximal, corresponding to max- 
imal mixing, when 1 1 = \j\[2. Though a careful mea- 
surement of P™ c may yield a value of |kp3| 2 slightly dif- 
ferent from 1/2, one here faces the difficulty of measuring 
a small term of order (1/^/2— |J7 M 3|) 2 . The vacuum term 
being quadratic in the deviation from maximality of muon 
neutrino flavor mixing contrasts strikingly with the mat- 
ter effect term being linear in the said deviation. 

Let us now consider JJJ), keeping terms only to lin- 
ear order in A. The A 2 terms will be shown to cancel 
from our effect and the corrections will be shown to be of 
0[{U e3 A) 3 }. Note that, for the MINOS baseline length 
and energies, |^4/ A31 1 = O(10 _1 ). It is useful to note 
that 



*l m =sl 3 (l + 2AA- 1 1 c 2 13 ) + (D(A 2 ), 
c 2 gm ^c 2 3 (l-2AA^ s 2 13 )+0(A 2 ), 

^^±^L = 0(A 2 ), 

. 2 BL . 2 A 31 L AL nn . , A rx 
sm — =sin — ^- cos 2#i3 sm(A 3:L L) 

2\ 



(9a) 
(9b) 

(9c) 



-0(A 2 ), 



sin 



. 2 A 3 i+A + P . 2 A31L AL 2 

-L = sm — ^— + — s 13 sin(A 3 iL) 

+0(A 2 ). 



(9d) 

L) 
(9e) 



Utilizing (9a-e) in JJJ), we obtain 
^ = + 24,^4,(4- 4, cos20 13 ) 
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FIG. 1: APjj^ as a function of the neutrino energy. The 
thick lines show the approximate analytic form given in Eq. 
IIUI I while the thin lines are obtained by the exact numerical 
solution of the equation of motion of the neutrinos traveling in 
matter. The comparison between our analytic results and the 
exact numerical solution is shown for three different values of 
| t/ e 3 1 2 . For the exact numerical solution we have used Sm^ = 
7 x 10 -5 eV 2 and sin 2 6 12 = 0.3. 
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FIG. 2: Difference in the number of neutrino and antineutrino 
events due to the deviation of \U^3\ from maximality. 
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Ea. l|l(J|l is the key result. It shows that the linear A 
term in AP W is proportional to 1 — 2|[/ At3 | 2 which is 
the deviation from maximality of the flavor mixing in 
vacuum, as evident from (JHJ). We may also note that the 
0(A 2 ) terms cancel between and Pp.fi- Moreover, 
the corrections involve 0[(Z7 e3 A) 3 ] terms since the low- 
est order A— terms in the transition amplitude [20] come 
with coefficients |f7 e3 | 2 and {7 e3 . Thus the A 3 terms are 
further suppressed. 

Our two significant approximations, i.e. ignor- 
ing ^m^-driven subdominant oscillations and the 
0[(C/ e3 A) 3 ] terms in Hl()(l . were for the convenience of 
analytical calculations. A numerical code has been de- 
veloped [2l| to calculate AP W , treating both the sub- 
dominant and the earth matter effect exactly. In Fig. ^ 
the thin (thick) lines show the AP MM of exact numerical 
evaluation (of 1|10|0 against E v with the other param- 
eters specified in the figure labels. The differences are 
quite small, thus enhancing our confidence in the use of 
((TUj) as a measure of |?7 e3 | 2 (l — 2|[/ /i3 | 2 ). Interestingly, 
even for |Z7 e 3| 2 = 0.01 when the magnitude of AP MM is 
a lot less, these differences are small. The explanation 
why subdominant oscillations are suppressed is that the 
AA21 term in the transition amplitude [20] comes with 
a coefficient proportional to £/ e3 . We have also checked 
numerically that, while the location (in E u ) of the AP W 
peak is sensitive to variations in Sm^, its magnitude is 
not. 

The sensitivity of MINOS to AP W is best discussed 
in terms of AiV^, the difference - due to the deviation 
of I 1 2 from maximality - between the number of ex- 
pected neutrino and antineutrino events. For an antici- 
pated MINOS exposure of 16 x 10 20 primary pro- 
tons on target (p.o.t), this has been plotted in Fig. 
along with the la statistical errorbars for the parameter 
values shown in the labels. Given the detection cross- 
section for D^s to be about half that of v^'s, we have 
assumed twice as much exposure for D^s as for v^s. The 
plot demonstrates the feasibility of such a measurement. 

In conclusion, we have shown how the measurement of 
AP Mfl , in a long baseline experiment such as MINOS, can 
probe the deviation l/\/2— lU^] from maximality of the 
flavor mixing of the muon neutrino. The exclusion (with 
errors taken into account) of a vanishing value of AP W 
will simultaneously demonstrate (7 e3 ^ ^ ^= — |L^, 3 |- 
The accuracy in the determination of |L7 e3 |(-l= — | C/ m3 |) 
will depend more sensitively on the measurement error 
in AP MAI than on the then uncertainty in the knowledge 
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of Sm^. 
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